Phase transition to the state with nonzero average helicity in dense neutron matter 
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The possibility of the appearance of the states with a nonzero average helicity in neutron matter 
is studied in the model with the Skyrme effective interaction. By providing the analysis of the self- 
consistent equations at zero temperature, it is shown that neutron matter with the Skyrme BSkl8 
effective force undergoes at high densities a phase transition to the state in which the degeneracy 
with respect to helicity of neutrons is spontaneously removed. 
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The issue of spontaneous appearance of spin polarized 
states in nuclear matter is a topic of a great current in- 
terest due to relevance in astrophysics. In particular, the 
scenarios of supernova explosion and cooling of neutron 
stars are essentially different, depending on whether nu- 
clear matter is spin polarized or not. On the one hand, 
the models with the Skyrme effective nucleon-nucleon 
(NN) interaction predict the occurrence of spontaneous 
spin instability in nuclear matter at densities in the range 
from go to 4gn for different parametrizations of the NN 
potential [H-Q (go — 0.16 fm -3 is the nuclear saturation 
density). On the other hand, for the models with the 
realistic NN interaction, no sign of spontaneous spin in- 
stability has been found so far at any isospin asymmetry 
up to densities well above go [To|]-[l6(. ^ n or d er to recon- 
cile two different approaches, based on the effective and 
realistic NN interactions, recently a new parametriza- 
tion of the Skyrme interaction, BSkl8, has been pro- 
posed fl7j , aimed to avoid the spin instability of nuclear 
matter at densities beyond the nuclear saturation den- 
sity. This is achieved by adding new density-dependent 
terms to the standard Skyrme force. The advantage of 
the BSkl8 parametrization is that it also preserves the 
high-quality fits to the mass data obtained with the con- 
ventional Skyrme force as well as it satisfactorily repro- 
duces the results of microscopic neutron matter calcula- 
tions (equation of state [lU, 1 5o pairing gap Hence, 
this Skyrme parametrization has a good potentiality in 
the studies of various neutron star properties [2(j ■ 

In terms of Landau Fermi liquid parameters, the fer- 
romagnetic instability of neutron matter is prevented 
by the requirement Go > —1, where Go is the zcroth 
coefficient in the expansion of the dimensionless spin- 
spin interaction amplitude on the Legendre polynomi- 
als, G = Y,i GiPi{cos9) Although this condition 
can hold true for all relevant densities, nevertheless, the 
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first coefficient Gi, with increasing density, can become 
large and negative, so that the condition Gi < — 3 could 
be reached in the high-density region of neutron mat- 
ter. These conditions, considered together, mean quite 
strong attractive interaction between neutron spins in the 
triplet spin state and repulsive, or weakly attractive in- 
teraction in the singlet spin state. As a result, neutron 
matter becomes unstable at high densities, and a new 
state characterized by a nonzero average helicity of neu- 
trons, (<xp°) ^ 0, is formed. Such a possibility was first 
studied with respect to an electron liquid in metals in 
Ref. [22j and, later, in the framework of a microscopic 
model, in Ref. [23| . Our primary goal here is to develop 
the proper formalism for the description of the states 
with a nonzero average helicity in neutron matter and to 
provide the corresponding analysis of the self-consistent 
equations for the BSkl8 Skyrme force. 

The nonsuperfluid states of neutron matter are de- 
scribed by the normal distribution function of neutrons 
/kik 2 = Tr pa+ 2 a Kl , where n = (p,a), p is momentum, 
a is the projection of spin on the third axis, and g is the 
density matrix of the system [H, Q . The self-consistent 
matrix equation for determining the distribution function 
/ follows from the minimum condition of the thermody- 
namic potential 0, [H[ and is 

/ = {exp(r e + F 4 ) + I}" 1 = {exp(YoO + I}" 1 ■ (1) 

Here the single particle energy e and quantity Y4 are 
the matrices in the space of k variables, with l4 KlK2 = 
Y4(5 KlK2 , Y"o = and Kj = — /.io/T being the Lagrange 
multipliers, /io being the chemical potential of neutrons, 
and T being the temperature. 

For the BSkl8 interaction, the ferromagnetic instabil- 
ity is avoided by adding to the standard Skyrme force 
the additional density-dependent terms. Nevertheless, 
although spontaneous spin polarization is excluded at all 
densities relevant for neutron stars, there is still the possi- 
bility, related to the spontaneous appearance of the state 
with a nonzero average helicity 
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A^<<rp°), p° = p/p, 



(2) 
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where (...) = tr /.../tr /, tr ... being the trace in the space 
of k variables. For this state, the single particle energy 
of neutrons reads 



e(p) = e (p)<j - A(p)erp°, 



(3) 



where <Ji are the Pauli matrices in the spin space and 
2A(p) is the energy splitting between the neutron spectra 
with different helicities (positive and negative). It follows 
from Eqs. Q) and <j3j) , that the distribution function of 
neutrons has the structure 



/(p)=/o(pVo + /||(pW°, 



(4) 



where 



fa = -{n(w+) + n(w_)}, /| 



■{n(u+) - n(u)-)}. (5) 



Here n{ui) = {exp(Y uj) + 1} 1 and 
w± = £o - Mo T A. 



(6) 



The quantity w± , being the exponent in the Fermi dis- 
tribution function n, entering Eqs. ([5]), plays the role 
of the quasiparticle spectrum. There are two branches 
of the quasiparticle spectrum, corresponding to neutrons 
with definite helicity, crp° = ±1. 

Note that the distribution function /o should satisfy 
the normalization condition 



y$Z/o(p) = 0. 
p 



(7) 



where g is the total density of neutron matter. The av- 
erage helicity A plays the role of the order parameter of 
a phase transition to the state, in which the majority of 
neutron spins are oriented along, or opposite to their mo- 
menta. By calculating the traces in Eq. @, it is easy to 
find that 



A 



E p /ii(p) 
E p /o(p)' 



(8) 



In order to get the self-consistent equations for the 
components of the single particle energy, one has to set 
the energy functional of the system, which reads 0, l24j 



E(f) = E (f) + E int {f), 
£ (/)=2^£ (p)/o(p), e (p) = 



(9) 



2m 



Eint(f) = ]T{£0(P)./0(P) + Ei(p)/i(p)}, 
P 



where 



eo(p) = ^E C/ o(k)/o(q), k= (10) 
q 



Here £ (p) ^ s the ^ ree sm gl e particle spectrum, too is 
the bare mass of a neutron, [/^(k), C/f(k) are the nor- 
mal Fermi liquid (FL) amplitudes, and £o,£i are the FL 
corrections to the free single particle spectrum. Using 
equation 5E = tr e(f)5f [21j, we get the self-consistent 
equations in the form 



£o(p) =£o(p) + £o(p) - Mo, 
&(p) = -A(p) Pi ° = ii(p). 



(12) 
(13) 



To obtain numerical results, we utilize the BSkl8 
parametrization of the Skyrme interaction, developed 
in Ref. 17j and generalizing the conventional Skyrme 
parametrizations. In the conventional case, the ampli- 
tude of NN interaction reads [2Q| 



v(p, q) = t (l + x a P a ) + -t 3 (l + x 3 P a )g a 
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(14) 



-^^(1 + Xl P a ){p 2 + q 2 



t 2 

^■(l + x 2 P a )pc[, 



where P a = (1 + <ti<T2)/2 is the spin exchange operator, 
ij, ajj and a arc some phcnomcnological parameters spec- 
ifying a given parametrization of the Skyrme interaction. 
The Skyrme interaction suggested in Ref. [ij) has the 
form 



w'(p,q) =«(P,<0 + ^2*4(1 + x 4 P CT )( P 2 + q 2 ) (15) 



fi 2 



i 5 (l + x 5 P a )pq. 



In Eq. (fT"5")) . two additional terms are the density- 
dependent generalizations of the t\ and t 2 terms of the 
usual form. 

The normal FL amplitudes U$,U\ can be expressed in 
terms of the Skyrme force parameters. For conventional 
Skyrme force parametrizations, their explicit expressions 
are given in Refs. [24|,[25j]. As follows from Eqs. (| X4[) and 
(|15l) , in order to obtain the corresponding expressions for 
the generalized Skyrme interaction (|15[) . one should use 
the substitutions 

t\ — > t\ + tigP, t\Xi — > t\Xi + t^x^gP , (16) 
t2 — > ti + hg 1 , t 2 x 2 -¥ t 2 x 2 + t 5 x 5 g' y . (17) 

Therefore, the FL amplitudes are related to the pa- 
rameters of the Skyrme interaction (|15[) by formulas 

C/ "(k) = 2< (1 - .to) + jg a (l - 13) + Jf[*i(1 - xi) 

(18) 

+ U(l - Xi )gP + 3t 2 (l + x 2 ) + 3*5(1 + T 5 )e 7 ]k 2 , 

[/f(k) = -2t (l - xo) - |e Q (l - x 3 ) + A[ t2 (i + X2 ) 

(19) 

+ i 6 (l + :r 5 )e 7 - ti(l - n) - t 4 (l - x 4 )/]k 2 
= o. + 6„k 2 . 
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It follows from Eqs. ([12]) and (JTSJl that 



where the effective neutron mass m n is defined by the 
formula 

7^— = ^— + ^Mi-xi)+U(i-x i )e fi (21) 

Zm n LrtiQ 8 

+ 3ta(l + S2) + 3*6(1 + 35)^], 

and the renormalizcd chemical potential [i should be de- 
termined from Eq. ([7]). 

Taking into account the explicit form of the FL ampli- 
tude U\ in Eq. (|T9)) . solution of Eq. (|T3l) for the energy 
gap A should be sought in the form 

A(p) = h -^up, (22) 

where v is some unknown quantity satisfying the equa- 
tion 

v = [Sn^ {q)dq - (23) 

This equation can be obtained from Eqs. (fTTj) , (|13[) af- 
ter passing from summation to integration, y Y] , . , — > 

I (2-1? h)* • ' •' an< ^ performing then the angle integration. 

Thus, with account of Eqs. ((5]) for the distribution 
functions /, we obtain the self-consistent equations (0 
and (|23j) for the renormalizcd chemical potential /i and 
the unknown is, determining the splitting A in the en- 
ergy spectrum ([3]) of neutrons with different helicities. 
Note that the energy spectrum ([3]) is invariant under the 
time reversion but not under the parity transformation. 
Hence, the state with A ^ is characterized by a spon- 
taneously broken P-symmetry. 

Now we present the solutions of the self-consistent 
equations at zero temperature for BSkl8 Skyrme 
force [l7|]- Note that the self-consistent equations have 
always the trivial solution A = (or v = 0), corre- 
sponding to the normal neutron matter. Besides, the 
self-consistent equations are invariant with respect to the 
change A — > —A, and, hence, nontrivial solutions for A 
enter in a pair with the same magnitude and opposite 
sign. The majority of neutrons will have positive helic- 
ity, if A > 0, and negative helicity, if A < 0. After 
solving the self-consistent equations, the average helic- 
ity A, playing the role of the order parameter of a phase 
transition, can be obtained from Eq. ([5]). According to 
Eq. both signs of the helicity of the given magni- 
tude are possible because of the two possible signs of the 
energy splitting A. Note that in order to preserve the 
realistic EoS of neutron matter obtained in Ref. [l8[ , the 
following constraints on the additional parameters of the 
BSkl8 parametrization were set 

= 7, t 4 (l-Xi) = -M 5 (l + x 5 ). (24) 
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FIG. 1. The energy splitting A(pf) between the neutron spec- 
tra with different helicities normalized to the neutron Fermi 
energy as a function of density at zero temperature for BSkl8 
interaction. 
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FIG. 2. Same as in Fig. [TJ but for the average helicity of 
neutron matter. 

Because of these constraints, the ti and is terms cancel 
completely in the FL amplitude Uo and in the effective 
neutron mass m n , and only the FL amplitude Ui is af- 
fected by the new terms. 

Fig. [U shows the energy splitting A(p = hkp) be- 
tween the branches of the neutron spectra with diffcr- 

h 2 k 2 

cnt helicities normalized to the Fermi energy £p = 2m F 
of the normal neutron matter as a function of den- 
sity. A spontaneous phase transition to the state with a 
nonzero helicity occurs at the critical density g w 5.67 go 
(go = 0.1586 fm~ 3 for BSkl8 force). The energy split- 
ting continuously increases with the density and becomes 
comparable with the neutron Fermi energy e p ■ Note that 
only the branch with the positive energy gap is shown in 
Fig. [1] while the symmetric branch (A — > —A) with the 
negative energy gap is not presented there. 

Fig. [5] shows the average helicity of neutron matter as 
a function of density obtained with the BSkl8 Skyrme 
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FIG. 3. Same as in Fig. [TJbut for the difference between the 
energy densities of the state with a nonzero average helicity 
and the normal state (A = 0) of neutron matter. 



interaction. The average helicity monotonously increases 
from zero till it is saturated and reaches the value A = 1 
at g w 6.93 go. Beginning from that density, all neutron 
spins will be aligned along their momenta (or opposite to 
them for the branch with the negative helicity, not shown 
in Fig.E]). 

In order to clarify whether the state with a nonzero 
average helicity is thcrmodynamically preferable over the 
normal state of neutron matter, we should compare the 
corresponding energies (at zero temperature). Fig. [3] 
shows the difference between the energy densities of the 
state with a nonzero average helicity and the normal state 
of neutron matter. It is seen that for all densities where 



nontrivial solutions (with A / 0) exist, this difference is 
negative and, hence, the state with the majority of neu- 
tron spins directed along (or opposite) to their momenta 
is preferable at that density range. 

In summary, we have considered the states with a spon- 
taneous nonzero average helicity in neutron matter with 
the BSkl8 Skyrme NN interaction, which are character- 
ized by broken parity. The self-consistent equations for 
the parameter, determining the energy splitting between 
the neutron spectra with different helicities, and the ef- 
fective chemical potential of neutrons have been obtained 
and analyzed at zero temperature. It has been shown 
that the self-consistent equations have solutions corre- 
sponding to a nonzero average helicity beginning from 
the critical density g ~ 5.67 go- Under increasing density, 
the magnitude of the average helicity increases and is sat- 
urated at the density g ~ 6.93 go, when all neutron spins 
arc aligned along (A = 1), or opposite (A = — 1) to their 
momenta. The comparison of the respective energies at 
zero temperature shows that the state with a nonzero av- 
erage helicity is preferable over the normal state at the 
densities beyond the critical one. The possible existence 
of the state with a nonzero average helicity in the dense 
core of a neutron star will affect the neutrino opacities, 
and, hence, may be of importance for the adequate de- 
scription of the thermal evolution of pulsars. In this re- 
spect, it is interesting to study also the impact of a strong 
magnetic field, characteristic for pulsars (magnetars), on 
the average helicity of dense neutron matter [22j ■ 
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